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Abstract 


Computational  fluid  dynamics  (CFD)  calculations  have  been  performed  for 
a  multi-body  system  consisting  of  a  main  missile  and  a  number  of 
submunitions.  Numerical  flow  field  computations  have  been  made  for 
various  orientations  and  locations  of  submunitions  using  an  unsteady,  zonal 
Navier-Stokes  code  and  the  chimera  composite  grid  discretization  technique 
at  transonic  speeds  and  zero  degree  angle  of  attack.  Both  steady  state  and 
unsteady  numerical  results  have  been  obtained  and  compared  for  two 
submunitions  and  a  missile  system.  Computed  results  show  the  details  of 
the  expected  flow  field  features,  including  the  shock  interactions. 
Computed  results  are  compared  with  limited  experimental  data  obtained  for 
the  same  configuration  and  conditions  and  are  generally  found  to  be  in  good 
agreement  with  the  data.  Comparison  of  the  unsteady  and  steady  state 
results  shows  an  appreciable  change  in  the  aerodynamic  forces  and 
moments. 
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COMPUTATIONAL  FLUID  DYNAMICS  MODELING  OF  MULTI-BODY 
MISSILE  AERODYNAMIC  INTERFERENCE 


1.  INTRODUCTION 

Aerodynamic  forces  and  moments  are  critical  design  parameters  used  in  the  design  of  artillery 
shells  and  bodies  flying  in  relative  motion  to  each  other.  The  advancement  of  computational  fluid 
dynamics  (CFD)  has  had  a  major  impact  on  projectile  design  and  development.1"4  Improved 
computer  technology  and  state-of-the-art  numerical  procedures  enable  scientists  to  develop 
solutions  to  complex,  three-dimensional  problems  associated  with  projectile  and  missile  aero¬ 
dynamics.  The  research  effort  has  focused  on  the  development  and  application  of  a  versatile 
overset  grid  numerical  technique  to  solve  multi-body  aerodynamic  problems.  This  numerical 
capability  has  been  used  successfully  to  determine  the  aerodynamics  on  a  multi-body  problem  of 
brilliant  antiarmor  (BAT)  submunition  dispersal  from  the  Army  tactical  missile  system  (TACMS). 
Figure  1  shows  a  schematic  diagram  of  this  multi-body  system. 


Figure  1.  Schematic  Diagram  of  Multi-bodv  System. 

The  complexity  and  uniqueness  of  this  type  of  multi-body  problem  result  from  the 
aerodynamic  interference  of  the  individual  components,  which  include  three-dimensional  (3-D) 
shock-shock  interactions,  shock-boundaiy  layer  interactions,  and  highly  viscous-dominated 
separated  flow  regions.  The  overset  grid  technique,  which  is  ideally  suited  to  this  problem,  involves 
generating  numerical  grids  about  each  body  component  and  then  oversetting  them  onto  a  base  grid  to 
form  the  complete  model.  With  this  composite  overset  grid  approach,  it  is  possible  to  determine 
the  3-D  interacting  flow  field  of  the  multi-body  system  and  the  associated  aerodynamic  forces  and 
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moments  at  different  positions  and  orientations  without  the  need  for  costly  regridding.  The 
solution  procedure  of  the  developed  technique  is  to  compute  the  interference  flow  field  at  multiple 
locations  until  final  converged  solutions  are  obtained  and  then  to  integrate  the  pressure  and  viscous 
forces  to  obtain  the  total  forces  and  moments.  The  complex  physics  and  fluid  dynamics  structure 
of  the  3-D  aerodynamic  interference  for  this  multi-body  problem  have  been  identified. 

A  description  of  the  computational  algorithm  and  the  chimera  technique  follows.  The  next 
section  describes  the  model  geometry  and  various  computational  grids  used  in  the  numerical 
computations.  Results  are  shown  for  both  steady  state  computations  for  the  missile  with  single  and 
multiple  BAT  submunitions  and  a  dynamic  computation  with  two  submunitions  at  transonic  speeds. 


2.  SOLUTION  TECHNIQUE 

2.1  Governing  Equations 

The  complete  set  of  3-D,  time-dependent,  generalized  geometry,  Reynolds-averaged,  thin 
layer  Navier-Stokes  equations  is  solved  numerically  to  obtain  a  solution  to  this  problem  and  can 
be  written  in  general  spatial  coordinates  T),  and  C,  as  follows:"5 

q  +  3^  F  +  3^  G  +  H  =  Re  S,  (1) 

in  which 

<;  =  Hfx,  y,  z,  t)  -  longitudinal  coordinate; 
r|  =  T|(x,  y,  z,  t)  -  circumferential  coordinate 
C,  =  C,(x,  y,  z,  t)  -  nearly  normal  coordinate; 
x  =  t  -  time 

In  Equation  1 ,  q  contains  the  dependent  variables:  density,  three  velocity  components,  and 
energy.  The  thin  layer  approximation  is  used  here,  and  the  viscous  terms  involving  velocity 
gradients  in  both  the  longitudinal  and  circumferential  directions  are  neglected.  The  viscous  terms 
are  retained  in  the  normal  direction,  and  are  collected  into  the  vector  S.  In  the  wake  or  the  base 
region,  similar  viscous  terms1  are  also  added  in  the  stream  wise  direction,  £.  An  implicit, 
approximately  factored  scheme  is  used  to  solve  these  equations. 

2.2  Numerical  Algorithm 

The  implicit,  approximately  factored  scheme  for  the  thin  layer  Navier-Stokes  equations  using 
central  differencing  in  the  rj  and  C,  directions  and  upwinding  in  £,  is  written  in  the  following  form:6 
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in  which  h  =  At  or  (At)/2  and  the  free-stream  base  solution  is  used.  Here,  5  is  typically  a  three- 
point  second  order  accurate  central  difference  operator,  8  is  a  midpoint  operator  used  with  the 
viscous  terms,  and  the  operators  8^  and  8^  are  backward  and  forward  three-point  difference 

operators.  The  flux  F  has  been  eigensplit,  and  the  matrices  A,  B,  C,  and  M  result  from  local 

linearization  of  the  fluxes  about  the  previous  time  level.  Here,  /  denotes  the  Jacobian  of  the 
coordinate  transformation.  Dissipation  operators  De  and  D.  are  used  in  the  central  space 

differencing  directions. 


2.3  Chimera  Scheme 

The  chimera  overset  grid  technique7-9,  which  is  ideally  suited  to  multi-body  problems, 
involves  generating  independent  grids  about  each  body  and  then  oversetting  them  onto  a  base  grid 
to  form  the  complete  model.  This  procedure  reduces  a  complex  multi-body  problem  into  a  number 
of  simpler  subproblems.  An  advantage  of  the  overset  grid  technique  is  that  it  allows  computational 
grids  to  be  obtained  for  each  body  component  separately  and  thus  makes  the  grid-generation 
process  easier.  Because  each  component  grid  is  generated  independently,  portions  of  one  grid  may 
lie  within  a  solid  boundary  contained  within  another  grid.  Such  points  lie  outside  the 
computational  domain  and  are  excluded  from  the  solution  process.  This  was  achieved  by 
modifying  Equation  2  for  chimera  overset  grids  by  the  introduction  of  the  flag  ib.  This  ib  array 
accommodates  the  possibility  of  having  arbitrary  holes  in  the  grid.  The  ib  array  is  defined  so  that  ib 
=  1  at  normal  grid  points  and  ib  =  0  at  hole  points.  Thus,  when  ib  =  1,  Equation  2  becomes  the 
standard  scheme,  but  when  ib  =  0,  the  algorithm  reduces  to  AQn  =  0  or  Qn+1  =  Qn  >  leaving  Q 
unchanged  at  hole  points.  The  set  of  grid  points  that  form  the  border  between  the  hole  points  and 
the  normal  field  points  are  called  inter-grid  boundary  points.  These  points  are  updated  by 
interpolating  the  solution  from  the  overset  grid  that  created  the  hole.  Values  of  the  ib  array  and  the 

interpolation  coefficients  needed  for  this  update  are  provided  by  a  separate  algorithm.7  Figure  2 
shows  an  example  where  the  parent  missile  grid  is  a  major  grid  and  the  BAT  submunition  grid  is  a 
minor  grid.  The  submunition  grid  is  completely  overlapped  by  the  missile  grid,  and  thus,  its  outer 


3 


boundary  can  obtain  information  by  interpolation  from  the  missile  grid.  Similar  data  transfer  or 
communication  is  needed  from  the  submunition  grid  to  the  missile  grid.  However,  a  natural  outer 
boundary  that  overlaps  the  submunition  grid  does  not  exist  for  the  missile  grid.  The  overset  grid 
technique  creates  an  artificial  boundary  or  a  hole  boundary  within  the  missile  grid  that  provides  the 
required  path  for  information  transfer  from  the  submunition  grid  to  the  missile  grid.  The  resulting 
hole  region  is  excluded  from  the  flow  field  solution  in  the  missile  grid. 


Missile  (Major)  Domain  Sixbmunition  (Minor)  Domain 


-  i'  iv 

vv..\  -vi/. 

jSb  '  •  1 

Artificial  Boundary  in  Parent  Domain 


Figure  2.  Inter-grid  Communications. 


The  chimera  method  depends  on  three  functions:  domain  comiectivity,  aerodynamics,  and 
body  dynamics.  The  aerodynamics  code  depends  on  the  domain  connectivity  code  to  supply 
hole  and  interpolation  information.  The  domain  comiectivity  code,  in  turn,  depends  on  the  body 
dynamics  code  to  supply  the  location  and  orientation  of  the  moving  bodies  relative  to  the  primary 
body.  Finally,  the  body  dynamics  code  depends  on  the  aerodynamics  code  to  provide  the 
aerodynamic  forces  and  moments  acting  on  the  moving  bodies.  For  moving  body  problems,  all 
grids  are  allowed  to  move  with  6  degrees  of  freedom  (DOF)  relative  to  an  inertial  reference  frame. 
Accordingly,  bodies  can  move  with  respect  to  others  without  the  necessity  of  generating  new  grids. 
With  this  composite  overset  grid  approach,  it  is  thus  possible  to  determine  the  aerodynamics 
associated  with  the  bodies  without  the  need  for  costly  regridding.  This  also  eliminates  potential 
accuracy  problems  attributable  to  severe  grid  stretching  used  by  many  other  techniques,  such  as  the 
zonal  blocked  grid  method  commonly  used  in  CFD. 
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3.  MULTI-BODY  PROBLEM  DESCRIPTION 


The  TACMS-BAT  multi-body  problem  involves  the  radial  dispensing  of  several  BAT 
submunitions  (see  Figure  1)  at  a  transonic  speed  and  thus  was  ideally  suited  for  the  numerical 
capability11  described  earlier.  The  3-D  radial  dispensing  of  these  submunitions  depends  on  the 
initial  ejection  velocity.  The  flow  field  is  complex  and  involves  3-D  shock-boundary  layer 
interactions,  and  TACMS-to-BAT  as  well  as  BAT-to-BAT  interactions.  Detailed  experimental 
or  theoretical  data  were  not  available  to  help  evaluate  the  submunition  dispensing  phenomenon 
for  the  entire  BAT  system,  and  thus,  the  numerical  solution  of  this  problem  was  initiated.  The 
chimera  solution  procedure  was  thus  used  to  determine  the  aerodynamic  interference  effects,  and 
CFD  was  brought  into  the  developmental  phase  of  the  BAT  program  to  ensure  successful 
dispensing  of  the  submunitions. 

The  missile  carries  13  submunitions;  the  first  10  outer  BAT  submunitions,  and  then  the 
three  inner  BAT  submunitions,  are  radially  dispensed  from  the  TACMS.  Once  released  from  the 
missile  bay,  the  self-guided  BAT  submunitions  autonomously  disperse  over  the  hostile  territory, 
use  their  sensors  to  detect  targets,  and  deliver  shaped  charged  warheads.  The  concern  here  is  the 
flight  dynamics  and  aerodynamics  of  the  dispensing  phenomenon.  Application  of  the  advanced 
CFD  modeling  technique  to  this  multi-body  dispensing  problem  was  to  provide  realistic 
simulation,  detailed  understanding  of  the  underlying  aerodynamic  interference  effects,  and  design 
information  that  can  lead  to  successfully  dispensing  the  BAT  submunitions  from  the  TACMS. 

4.  MODEL  GEOMETRY  AND  GRIDS 

An  advantage  of  the  chimera  technique  is  that  it  allows  computational  grids  to  be  obtained 
for  each  body  component  separately  and  thus  makes  the  grid  generation  process  easier.  Figure  3 
shows  a  computational  grid  for  the  complete  model,  including  the  missile  and  the  BAT 
submunitions.  Also  shown  here  are  the  sections  of  the  three-dimensional  BAT  computational 
grids  overset  onto  the  missile  grid.  Figure  4  shows  a  computational  grid  for  one  BAT 
submunition.  As  part  of  the  chimera  procedure,  this  BAT  grid  is  partially  cut  by  the  missile 
body  itself.  Similarly,  the  presence  of  the  BAT  submunition  cuts  a  hole  in  the  missile  grid  (see 
Figure  5).  The  missile  grid  consists  of  three  zones:  one  on  the  nose  region  ahead  of  the  cavity, 
one  in  the  cavity  itself,  and  the  third  one  aft  of  the  cavity  region.  Each  of  these  three  zones  is  a 
rectangular  grid. 

The  grid  around  the  submunition  consists  of  two  zones  (one  for  the  body  and  one  for  the 
base  region)  and  was  obtained  using  an  C-topology  and  a  rectangular  topology,  respectively.  The 
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submunition  grids  were  individually  generated  and  then  overset  as  shown  in  Figure  3  to  form  the 
complete  grid  system.  The  computational  grids  shown  here  correspond  to  the  pitch  plane.  The 
missile  grid  serves  as  the  main  background  grid  for  the  computations.  Figure  4  shows  a 
computational  grid  for  computations  with  the  BAT  submunition  at  a  distance  about  a  diameter 
away  from  the  center  line  of  symmetry  of  the  missile.  For  steady  state  or  unsteady  dynamic 
computations,  the  same  submunition  grids  are  used,  and  there  was  no  need  to  regenerate  new 
submunition  grids.  Typically,  the  missile  grid  consisted  of  983,000  points.  The  entire  grid 
system  consisted  of  about  1 .2  million  points  for  the  missile  and  the  one-B  AT  (with  and  without 
sting)  case  as  well  as  for  the  two-BAT  cases.  Note  that  the  grid  setup  allows  computation  of  the 
base  region  flow  field  of  the  submunitions.  Grid  points  are  clustered  near  the  missile  and  the 
BAT  submunition  surfaces  to  capture  the  viscous  boundary  layers.  No  attempt  has  been  made 
to  adapt  the  computational  grids  to  gradients  in  the  flow  field  variables. 


Figure  3.  Grids  for  the  BAT  Submunition  Dispensing  From  TACMS. 
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The  actual  cavity  surface  of  the  missile  bay  where  the  BATs  are  stored  in  their  original  positions 
was  used  in  the  CFD  computations.  The  flow  field  in  the  bay  is  viscously  dominated,  turbulent,  and 
quite  complex,  with  the  BAT  located  in  the  near  field.  It  is  difficult  to  accurately  determine  the 
interference  effects  by  theoretical  or  experimental  means.  This  is  especially  true  when  the  BAT  is 
submerged  in  the  bay.  Limited  wind  tunnel  experimental  data12  are  available  for  a  reduced  scale  model 
for  the  missile  and  BAT  submunitions.  However,  such  data  can  suffer  from  sting  effects  and  for 
viscously  dominated  cavity  flows,  they  may  not  scale  to  the  real  flight  conditions. 

5.  RESULTS 

Both  steady  state  and  a  dynamic  unsteady  numerical  calculations  have  been  performed  to 
numerically  simulate  the  missile  and  the  BAT  system.  Computations  have  been  run  at  =  1.2 
and  1.5.  Computational  modeling  is  restricted  to  the  symmetrical  submunition  dispersal.  Here,  the 
missile  is  at  0°  angle  of  attack,  and  the  BATs  are  dispensed  symmetrically  following  the  same  radial 
trajectory  away  from  the  projectile.  Appropriate  symmetry  is  used  for  the  one-BAT  and  multiple 
BAT  cases.  For  the  multiple  BAT  cases,  the  computational  domain  consists  of  a  36°  segment  in 
the  circumferential  plane  (see  Figure  6).  Also  shown  is  the  submunition  grid,  which  is  entirely 
contained  in  the  background  projectile  grid.  Because  of  symmetry,  the  requirements  for  grid  sizes, 
computer  resources  such  as  computer  memory,  and  run  time  are  reduced. 


Figure  6.  Circumferential  Cross-Sectional  Grid. 
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The  CFD  modeling  was  initially  used  to  compute  transonic  flow  over  the  wind  tunnel  model 
consisting  of  the  missile  and  one  sting-mounted  BAT.  Computations  were  also  performed  for  the 
BAT  without  the  sting  for  direct  comparison  (see  Figure  7).  Computed  Mach  contours  show  some 
influence  of  the  sting  on  the  solution  near  the  aft  end  of  the  cavity.  Quantitatively,  the  presence  of 
a  sting  affects  the  forces  and  moments  acting  on  the  BAT  submunition  and  to  a  lesser  extent,  on  the 
missile  itself.  Computed  surface  pressures  obtained  from  the  flow  field  solution  for  the  BAT 
submunition  were  compared  with  experimentally  measured  pressure  on  the  BAT  and  were  found 
to  be  in  very  good  agreement  for  two  different  radial  locations  of  the  submunitions  (see  Figure  8). 
Computed  surface  pressure  on  the  missile  cavity  surface  is  shown  in  Figure  9  for  the  one-BAT 
case.  Computed  surface  pressures  again  agree  well  with  the  experimental  data.  These  comparisons 
inspire  confidence  in  the  quality  of  the  CFD  results  for  the  multi-body  system  consisting  of  the 
missile  and  one  BAT.  The  CFD  modeling  includes  the  wake  and  the  base  region  of  the  BAT  as 
they  would  be  in  real  flight  situations. 

For  the  actual  system,  which  consists  of  the  TACMS  and  the  10  outer  BAT  submunitions, 
the  present  CFD  capability  was  then  used  to  compute  the  flow  field  for  the  flight  conditions. 

This  case  now  included  not  only  BAT-to-TACMS  interactions  but  also  3-D  BAT-to-BAT 
interactions.  Figure  10  shows  the  results  of  such  a  numerical  simulation  for  the  symmetrical 
dispensation  of  all  10  BAT  submunitions  from  the  TACMS.  The  computed  Mach  contours  are 
shown  at  two  selected  positions  along  the  BAT  submunitions.  The  3-D  BAT-to-BAT 
interactions  are  clearly  depicted  in  this  figure.  To  promote  better  understanding  of  the 
aerodynamic  interference  associated  with  the  TACMS-BAT  system,  numerical  computations 
were  performed  with  the  BAT  submunitions  at  various  radial  positions  away  from  the  missile. 
Figure  1 1  shows  the  computed  results  for  these  cases  for  the  actual  transonic  flight  conditions 
from  bottom  to  top.  The  first  position  corresponds  to  the  BAT  submunitions  semi-chambered 
in  the  bay  (cavity)  of  the  missile.  The  next  two  positions  are  approximately  20.4  and  25.8  inches 
away  from  the  centerline  of  the  missile.  These  pictures  clearly  show  the  complicated 
aerodynamic  interference  between  the  missile  and  the  BAT  submunitions  in  a  longitudinal  view. 

The  flow  inside  the  cavity  is  primarily  subsonic  and  contains  large  regions  of  separated  flow. 

The  dramatic  change  in  the  interacting  flow  field  can  be  observed  as  the  BAT  submunitions  move 
farther  and  farther  radially  out  from  the  semi-chambered  position.  The  flow  field  over  the 
submunitions  is  3-D,  and  the  computed  results  clearly  show  this  feature  in  the  wake  or  base 
region  of  the  submunitions. 
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Computations  were  also  performed  for  this  multi-body  problem  with  five  BATs  at  Mx  = 

1 .2  and  a  =  0°.  This  case  included  modeling  of  half  of  the  actual  missile  bay.  Because  of 
symmetry,  CFD  modeling  uses  four  side  BATs  and  a  half  BAT  on  the  top  (lee  side)  as  well  as 
the  bottom  (wind  side).  Figure  12  shows  the  Mach  number  contours  for  this  case.  Qualitatively, 
it  shows  the  expected  shock  structure  and  the  flow  field  resulting  from  the  submunition 
interactions.  Figures  13  and  14  show  the  circumferential  Mach  number  contours  for  the  five- 
BAT  case  at  two  longitudinal  stations,  2.5  and  3.2  calibers  from  the  nose  of  the  missile.  Both 
locations  are  in  the  cavity  of  the  missile.  These  figures  show  the  BAT-to-BAT  interactions  and 
the  effect  of  the  cavity  shape  on  the  solutions.  Figure  13  indicates  a  smaller  region  of  low  speed 
flow  between  the  cavity  surface  and  the  bottom  surfaces  of  the  BATs  and  high  speed  flow  near 
the  top  surfaces  of  the  BATs.  At  a  station  down  stream  in  the  cavity,  Figure  14  shows  a  large 
region  of  the  low  speed  flow  between  the  BATs  and  missile  cavity  surface  as  well  as  near  and 
away  from  the  top  surfaces  of  the  BATs.  Aerodynamic  forces  and  moments  were  obtained  from 
the  computed  solutions.  Figure  1 5  shows  normal  force  coefficient,  axial  force  coefficient,  and 
pitching  moment  for  the  submunitions  as  a  function  of  radius  (measured  from  the  centerline  of 
the  missile).  These  computed  force  (CN,  CA)  and  moment  (CM)  coefficients  were  compared 
with  the  experimental  data  and  are  found  to  be  in  good  agreement  with  the  data.  The  measured 
axial  force  coefficient,  CA,  does  not  include  base  drag  and  computed  axial  force  coefficient  also 
excludes  base  drag  of  the  submunitions  for  direct  comparison. 


Figure  12.  Mach  Contours  for  the  Five-BAT  Case. 


14 


Figure  14.  Mach  Contours  X/D  =  3.2  Calibers. 
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Figure  15.  Force  and  Moment  Coefficients. 


A  dynamic  simulation  was  performed  for  the  missile  and  the  two-BAT  case  at  =1.5 
and  a  =  0°.  The  dynamic  computations  were  started  from  a  converged  steady  state  solution  with 
the  submunitions  located  in  the  cavity  chamber.  The  same  computational  grids  (see  Figure  6) 
were  used  for  the  entire  dynamic  simulation  event  without  the  need  of  regridding.  The  chimera 
procedure  described  earlier  was  also  used  for  this  computation.  The  forces  and  moments 
provided  by  the  aerodynamics  code  were  used  by  the  body  dynamics  code  to  determine  the 
location  and  orientation  of  the  moving  submunitions  relative  to  the  missile.  This  required  domain 
connectivity  (hole  and  interpolation  information)  at  each  time  step.  The  domain  connectivity 
information  was  then  used  in  the  aerodynamics  code  to  provide  the  aerodynamic  forces  and 
moments.  This  procedure  was  repeated  at  each  time  step  during  the  dynamic  simulation.  Figure 
16  shows  three  snapshots  of  computed  pressure  contours  in  time.  It  shows  the  flow  field 
changing  significantly  with  the  submunition  being  dispensed.  The  bottom  picture  corresponds  to 
the  initial  position  of  the  submunitions.  The  other  two  correspond  to  16  and  45  msec  in  time. 

The  submunitions  were  restricted  to  the  pitch  plane  and  were  not  allowed  to  yaw  or  roll. 
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Figure  1 7  shows  the  force  and  moment  history  as  a  function  of  time.  Shown  here  are  the 
normal  force,  CN,  axial  force,  CA,  and  pitching  moment,  CM,  coefficients  for  the  submunition. 
Since  the  computations  include  two  BAT  submunitions,  BAT-to-BAT  interactions  are  included. 
These  interactions  are  critical  and  have  a  strong  effect  on  the  aerodynamic  forces  and  moments. 
The  normal  force  and  pitching  moment  coefficients  indicate  the  unsteady  nature  of  the  interacting 
flow  field.  A  case  that  corresponds  to  16  msec  in  time  during  the  dynamic  simulation  was 
selected  and  frozen  to  run  a  steady  state  computation.  This  case  corresponded  to  the 
submunition  being  at  an  angle  -1°  nose  down  (pitch  angle).  A  static  converged  solution  was  then 
obtained  for  the  missile  and  the  submunitions  in  that  orientation.  Comparison  of  the  static  and 
dynamic  results  for  this  case  is  shown  in  Table  1 .  Again,  all  three  force  and  moment  coefficients 
are  included.  The  axial  force  coefficient  (with  or  without  base  drag)  does  not  change  appreciably; 
however,  normal  force  and  pitching  moment  coefficients  are  seen  to  change  significantly.  These 
results  indicate  the  importance  of  the  dynamic  numerical  simulation,  which  may  be  needed  to 
accurately  predict  the  aerodynamics  of  the  submunition  dispersal. 


Figure  17.  Time  History  of  Force  and  Moment. 
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Table  1 


Static  and  Dynamic  Force  and  Moment  Comparison 


DYNAMIC 

STATIC 

R  =  20.23  fsi 

R  =  20.23  fsi 

CN 

+  0.187 

+0.393 

Cm 

- 1.026 

-  0.790 

CA 

w/base  drag 
w/o  base  drag 

+  0.776 
+  0.560 

+  0.794 
+  0.586 

6.  CONCLUDING  REMARKS 

A  computational  study  was  undertaken  to  compute  the  3-D  flow  fields  for  a  multi-body 
system  consisting  of  a  missile  and  multiple  BAT  submunitions.  Flow  computations  were 
performed  at  transonic  speeds  =  1.2, 1 .5)  and  a  =  0.0°  using  a  3-D  unsteady  Navier-Stokes 
code  and  chimera  composite  grid  discretization  technique.  Overset  body-conforming  grids  were 
used  to  individually  model  the  missile  and  the  BAT  submunitions.  Computed  results  have  been 
obtained  for  different  locations  of  the  submunitions.  Computed  pressure  and  Mach  contours 
show  the  details  of  the  3-D  aerodynamic  interference  flow  field  for  the  missile  and  the 
submunitions.  The  computed  flow  field  includes  both  the  missile-to-BAT  as  well  as  BAT-to- 
BAT  interactions.  Both  computed  results  for  sting-mounted  and  no-sting  models  have  been 
obtained  for  the  one-BAT  case,  and  the  computed  results  do  not  show  appreciable  change  in  the 
surface  pressures.  Computed  results  for  the  multiple  BAT  cases  include  BAT-to-BAT 
interactions  and  show  their  large  effect  on  the  computed  forces  and  moments.  Computed  surface 
pressures  as  well  as  forces  and  moments  have  been  compared  with  the  experimental  results  for 
the  same  configuration  and  conditions  and  are  generally  in  good  agreement  with  the  data. 
Dynamic  unsteady  results  were  also  obtained  for  a  multiple  BAT  case  and  compared  with  the 
static  results  for  the  same  location  and  orientation  of  the  submunitions.  These  comparisons 
show,  in  some  cases,  large  changes  in  the  aerodynamic  force  and  moments  between  the  static  and 
dynamic  computations.  Future  study  will  include  modeling  of  asymmetrical  submunition 
dispersal,  which  will  require  full  3-D  computations  and  large  computing  resources. 
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This  work  represents  the  application  of  a  chimera  overlapping  grids  approach  for  accurate 
numerical  calculation  of  aerodynamics  involving  multiple  bodies  with  and  without  relative 
motion.  The  predictive  numerical  capability  has  been  used  to  provide  the  development 
community  numerical  data  and  basic  flow  field  design  information  to  more  effectively  guide  the 
design  of  a  multi-body  missile  configuration.  It  allows  accurate  and  realistic  numerical  prediction 
of  interference  effects  and  aerodynamics  required  for  the  improved  design  and  modification  of 
current  and  future  multi-body  missile  and  projectile  configurations. 
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speeds  and  zero  degree  angle  of  attack.  Both  steady  state  and  unsteady  numerical  results  have  been  obtained  and  compared  for 
two  submunitions  and  a  missile  system.  Computed  results  show  the  details  of  the  expected  flow  field  features,  including  the 
shock  interactions.  Computed  results  are  compared  with  limited  experimental  data  obtained  for  the  same  configuration  and 
conditions  and  are  generally  found  to  be  in  good  agreement  with  the  data.  Comparison  of  the  unsteady  and  steady  state  results 
shows  an  appreciable  change  in  the  aerodynamic  forces  and  moments. 
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